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Abstract 

Two-dimensional spatially localized structures in the complex Ginzburg-Landau equation with 1:1 
resonance are studied near the simultaneous occurrence of a steady front between two spatially 
homogeneous equilibria and a supercritical Turing bifurcation on one of them. The bifurcation 
structures of steady circular fronts and localized target patterns are computed in the Turing-stable 
and Turing-unstable regimes. In particular, localized target patterns grow along the solution 
branch via ring insertion at the core in a process reminiscent of defect-mediated snaking in one 
spatial dimension. Axisymmetric oscillons on these solution branches are found to be stable over 
a wide parameter interval, and subject to various types of instability otherwise. Direct numerical 
simulations reveal novel depinning dynamics of localized target patterns in the radial direction, and 
of circular and planar localized hexagonal patterns in the fully two-dimensional system. 
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1. Introduction 

For several decades, two-dimensional (2D) localized structures have been a frequent subject of 
experimental, numerical and analytical studies [29|. A particularly interesting class of localized 
structures, known as localized states (LS), results from the coexistence of two stationary states, 
either of which can fill the entire spatial domain [l6|. These states may take the form of spatially 
homogeneous equilibria or spatially periodic patterns. A stable pattern can be created from an 
equilibrium Aq by the Turing instability, which can be either supercritical or subcritical. In the 
supercritical case, since the equilibrium Aq is unstable while the pattern is stable, the latter invades 
the former via a pulled front and thus no steady LS are expected [3l| ; Ref. 0 describes the invasion of 
a homogeneous state by a hexagonal pattern of this type. In the subcritical case, Aq and the pattern 
can both be stable and the front between them (known as a pushed front) may remain pinned to the 
pattern over a range of parameter values [28|, • This pinning mechanism underlies much recent 

work on both ID and 2D LS, notably in the quadratic-cubic Swift-Hohenberg equation (SH2^. 
The bifurcation structures and temporal dynamics of ID LS in SH23 are now well understood [3[: 
stationary reflection-symmetric LS are located on a pair of intertwined branches, one of which 
corresponds to states with maxima in the center and the other to states with minima in the center. 
As one proceeds along either branch these states grow in length via the repeated nucleation of 
a pair of wavelengths, one on either side of the LS, in a process reflected in the back-and-forth 
oscillations of the branches and referred to as standard homoclinic snaking. The two LS branches 
are in addition interconnected by an infinite number of cross-links (known as “ladders”) formed of 
stationary but asymmetric LS. These states, together with the associated multipulse states, coexist 
within a parameter interval referred to as the snaking or pinning region. Outside this interval the 
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LS grow or decay in time via a process referred to as depinning, i.e., through the nucleation or 
destruction of cells at the location of the fronts in a time-dependent manner, while leaving the 
interior cells undisturbed. See 0 for a recent review. 

More recently, 2D LS in SH23 have also been investigated. In Ref. [ 2 ^ localized d-dimensional 
axisymmetric solutions, hereafter referred to as localized target patterns, are studied using numerical 
continuation. The time-independent axisymmetric SH23 in d dimensions is treated as a spatial 
dynamical system in the radial coordinate r, with d treated as an additional continuous parameter. 
For localized target patterns with a large radius, the radial profile of the front between the target 
and the equilibrium Aq resembles a ID pushed front. However, in contrast to ID, near r = 0 (known 
as the core region) the target typically has a larger amplitude than at large r (known as the far 
field). As a result new rings can be nucleated both in the core and at the front, implying that the 
snaking behavior of localized target patterns for d > 1 will differ from standard homoclinic snaking 
at d = 1. In particular, when the radius is sufficiently large, new rings are nucleated exclusively in 
the core. A different type of localized structure, stationary spatially localized hexagonal patterns, 
are studied in Ref. [HI using numerical continuation. The existence of these solutions follows from 
a spatial conservation law for the time-independent SH23 in 2D that generalizes a conserved spatial 
Hamiltonian present in ID. Two classes of such states are known, hereafter referred to as planar 
and fully localized hexagonal patterns. The former are localized in ID, forming a hexagon-filled 
stripe connected to the background state by a pair of planar fronts, while the latter are bounded by 
a hexagonal front and localized in 2D. Both structures grow via the nucleation of additional cells 
along the bounding fronts, in a process similar to standard homoclinic snaking. A novel feature of 
the planar case is the presence of secondary snaking branches, corresponding to the nucleation of 
individual cells along either front, between successive oscillations of the primary snaking branch that 
reflect the addition of complete rows of cells along the fronts. This growth mechanism is also observed 
along the snaking curve of fully localized hexagons. Outside the snaking region time-dependent 
calculations show that these almost planar solutions evolve by losing or gaining cells one-by-one in 
a manner that recapitulates their growth along the snaking branch of time-independent solutions. 

In the above example one of the two states comprising the LS is homogeneous in space while the 
other is a steady, spatially periodic Turing pattern. In other PDFs the two competing states may be 
both homogeneous states or both stationary patterns. A notable pattern-forming PDE exhibiting 
bistability between homogeneous states is the forced complex Ginzburg-Landau equation (FCGLE) 
with 2:1 resonance. This is an equation for a complex scalar amplitude A = U-\-iV [5| which describes 
the effect of subharmonic forcing on a Hopf bifurcation to a homogeneous oscillation in spatially 
extended systems [l|. This PDE possesses an up-down symmetry A —A and exhibits a parameter 
regime where the two “equivalent” equilibria TAq related by this symmetry are simultaneously 
stable. In ID a codimension-zero family of steady Ising fronts between TAq is generically present, and 
may undergo a pitchfork bifurcation into traveling Bloch fronts Q ; near this Ising-Bloch transition 
the domain walls can exhibit rich spatiotemporal dynamics In Ref. 0, a dynamical theory 
for the domain wall between ±Ao in 2D leads to a characterization of different growth regimes and 
the prediction of the existence of stable “droplets”. In this paper we study instead the 1:1 EGGLE 

BA 

-^ = (1 + ia)V‘^A + (/i + in)A - (1 + il3)\A\^A 7, (1.1) 

which exhibits, under appropriate conditions, bistability between two equilibria A+ = 
and A~ = U~ ^iV~ belonging to the upper and lower branches of an S-shaped bifurcation diagram 
(Eig. [ 1 ]). In ID this regime generically exhibits a codimension-one family of steady Bloch fronts 
between A^ when both A^ are stable. Steady LS known as type-I LS can be constructed by 
assembling two such fronts back-to-back. Recent work has focused on new dynamical phenomena 
that arise when A+ becomes Turing-unstable, including the sna king behavior and depinning 
dynamics of steady localized Turing patterns known as type-H LS [23|. In contrast to standard 
homoclinic snaking, type-H LS grow along a single snaking branch by nucleating new cells at the 
center of the wavetrain in a process termed defect-mediated snaking (DMS) Outside the snaking 
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Figure 1: Sketch of the bifurcation diagrams of the building blocks of (a) ID and (b) 2D LS in the 1:1 FCGLE (11.11) . 
Two equilibria, labeled and separated by a pair of saddle-node bifurcations 7 ^“^, may be simultaneously stable. 
The Turing bifurcation on A+ at 7 ^ produces a branch of stripes S in ID, and branches of stripes S and hexagons 
H in 2D. The solid (dashed) segments are stable (unstable). 


region, type-II LS depin by creating or destroying cells in the interior of the wavetrain via successive 
phase slips, as studied quantitatively in Ref. [25|. In the present paper, we characterize 2 D LS in 
the (x^y) plane in the parameter regime where type-I and type-II LS are found in ID. 

The bifurcation diagrams of the building blocks of these ID and 2 D LS are shown schematically 
in Fig.[T]in terms of a suitable norm \\A\\ as a function of the control parameter 7 . In ID, a branch of 
stripes S bifurcates supercritically from ^4+ at the Turing bifurcation 7 ^, and may become bistable 
with A~ to form type-II LS (Fig. [Ha)). In 2 D, the Turing bifurcation generates, in addition to 
S', a pair of hexagonal patterns, hereafter referred to as that differ in whether the hexagons 
superposed on ^4+ are composed of bumps or holes [ll| . As a result the stripe pattern is unstable 
to hexagonal perturbations even when it bifurcates supercritically. The bifurcate transcritically 

and are likewise both unstable. However, in general one or the other of the folds back at 7 ^ > 7 ^ 
and acquires stability as indicated in Fig. [Hb). This is the case in the present case as well. 

We begin by focusing on one-pulse solutions in which the region occupied by the A+ state forms a 
single connected component embedded in the A~ state, although we briefly comment on “reciprocal” 
localized states in which the A~ state is embedded in the A+ state. When the connected component 
is finite, the interface separating A^ forms a closed curve, and we refer to the region inside (outside) 
this curve as the inner (outer) region. For most of the localized states considered in this paper, 
this curve takes the shape of a perfect circle, the only shape that respects the rotational symmetry 
of Eq. (|l.ip . This behavior is reminiscent of variational systems that tend to minimize a surface 
energy, even though Eq. (HID is nonvariational. In ^we study the simplest case, the case of circular 
fronts between A^. Next, we consider cases where the A+ state is Turing-unstable, leading to a 
patterned state near A+. As shown in 0the steady-state solutions that result may take the form of 
an (approximately) periodic structure in the radial direction, and we find localized target patterns 
of this type using numerical continuation. In ^4.11 we address the more complicated case where the 
pattern has hexagonal structure, and use direct numerical simulation (DNS) to identify a new type 
of LS consisting of localized hexagons that are surrounded by a ring, with this compound structure 
embedded in A~. We refer to these states as circular localized hexagonal patterns. 

In order to study axisymmetric solutions, we write the 1:1 ECGLE (HU in polar coordinates 

(r, 0), 

At = {1 + ia) (drr H—7 H— A + f{A), (1-2) 


where f{A) = (/i -h in)A — (1 + il3)\A\‘^A -t- 7 . Steady axisymmetric solutions satisfy Eq. (jl. 2 p with 
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dtA = deA = 0 , i.e. 


(1.3) 


(drr + -dr) A + = 0- 

Equation can be readily generalized from 2D to d dimensions, 

(drr + ^dr)A+r^f{A)=0, (1.4) 

\ r J la 

where d is treated as a continuous parameter (iil. The boundary conditions for axisymmetric 
solutions are taken to be Neumann boundary conditions (NBC), i.e., = 0 at r = 0 , oo, 

independently of d, or in terms of the real and imaginary parts of A, 

Ur = Vr = 0, r = 0,oo. (1.5) 

The L^-norm of a d-dimensional axisymmetric solution A{r) on a circular domain r G [0, R] is 
measured using either the ID radial norm 


Piu 



/ 


\A{r)\‘^dr^ 


( 1 . 6 ) 


or the d-dimensional radial norm 


r,d • — 





idr 



rd 


(1.7) 


Although the d-dimensional radial norm is the physical L^-norm for a multi-dimensional solution, 
the ID radial norm weights all parts of the solution profile equally and hence better represents the 
profile changes near r = 0. The L^-norm for a 2D solution A(x, y) (not necessarily axisymmetric) 
on a 2D domain Q is defined as 


Pll(a;,y) := JJjA{x,y)\^dxdy /li dxdy. (1-8) 

This fully 2D norm is compatible with the d-dimensional radial norm because Eq. dm with d = 2 
is identical to Eq. when A depends on r only and ^2 is a circular domain of radius R. Eor a 
localized state consisting of an inner pattern A* of radius p embedded in an outer pattern A^ on a 
domain of radius R {R> p)^ any of the above L^-norms generates a monotonic function of p that is 
well approximated by 

= + (1-9) 

The spectral stability of an axisymmetric LS AQ{r) depends on its stability with respect to 
perturbations proportional to or equivalently with azimuthal wavenumber m, where m G Z. 

Hereafter we refer to this perturbation via its symmetry group as a perturbation, where B^ 
denotes the dihedral group of order m. Thus Bq represents continuous rotations, Bi represents the 
trivial rotation by 27r, while B^ for m > 2 is generated by rotations by 27r/m together with the 
reflections 0 —0. The B^ stability of Ao(r) is given by the spectrum of the linearized ECGL 

operator with dee = namely 

3/72 + Vg 2U0V0 

2U0V0 + 3 Eo' J ’ 

with the boundary conditions given by Eq. m- 
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In addition to the above fully localized states in which either the state or the patterned state 
occupies a finite region, there exist planar localized states in which these states occupy a region that 
is infinite along a certain direction. If such localized states are computed on a rectangle with periodic 
boundary conditions (PBC), we may suppose that the solutions are localized in, say, the x-direction, 
and periodic in the orthogonal direction, namely the ^-direction. Two solutions of this type, planar 
fronts between ^4+ and A~ and localized stripe patterns can be constructed by extending Type-I LS 
and Type-II LS in the ^-direction. Because these states are invariant under translation in both 
satisfy the same PDE as their ID counterparts and so will not be discussed further. Instead we 
focus in g4.3l on nontrivial examples in which the inner pattern is itself periodic in 2D. We conclude 
in m with some open questions. 


As in Refs. [2J, [23 we fix three of the parameters in Eq. II.li viz., a = —1.5, /3 = 6 and // = —!, 


and vary the remaining two parameters u and 7. Numerical continuation and DNS are, respectively, 
performed using AUTO-O7p @ and the ETD2 scheme Unless otherwise specified, the domain 
size for DNS is [—100,100] x [—100,100]. 


2. Circular fronts 


The properties of ID fronts between in Eq. (HH) are well understood in Hi following 
earlier work in [l8|. In a bifurcation diagram of Type-I LS formed by two steady ID fronts 
assembled back-to-back, the successive folds converge on a limiting parameter value in an oscillatory 
fashion as A+ grows longer [23 • This bifurcation behavior is termed collapsed snaking^ and the 
limiting parameter value is termed the collapse point. At fixed v = b the collapse point occurs at 
7^*^ = 1.8419. Traveling ID fronts at 7 7^ 7^*^ correspond to heteroclinic orbits between A^ of an 
ODE in the comoving frame (iBj . 

In this section we study a family of 2D circular fronts between A^ in Eq. (HH). As in Refs. Hi HI 
we fix z/ = 5 and vary 7 as the bifurcation parameter. In ^2.1l we compute the bifurcation diagram 
of steady circular fronts, and compare the observed bifurcation behavior with ID collapsed snaking. 
In ^2.21 we compute the spectral stability of this solution branch, and demonstrate the existence of 
stable oscillons and the presence of different types of symmetry-breaking instabilities. In ^2.31 we 
compute traveling circular fronts using both 2D and axisymmetric time evolution, and show that 
these correspond to heteroclinic orbits between A^ in the radially comoving ODE. 


2.1. Radially collapsed snaking 

As pointed out in ^ steady d-dimensional axisymmetric solutions to Eq. (HU satisfy Eq. diu. 
Since in ID (d = 1) the collapsed snaking branch of Type-I LS originates from the lower saddle-node 
^SN S-shaped bifurcation diagram of equilibria (Eig. [T]), we may postulate the same for the 

branch of steady circular fronts in 2D (d = 2). In ID weakly nonlinear solutions near 7;?^ can be 
derived using standard asymptotic methods (eg.. Ref. Q) and these can be continued from ID to 
2D by incrementing d from d = 1 to d = 2 in Eq. (lU; see Ref. for an analytical description of 
these 2D weakly nonlinear solutions in the 2:1 ECGLE. Once d = 2 has been reached the resulting 


solutions can then be continued in 7 with d = 2 fixed, namely within Eq. (HSJ. Eigure |2(a)| shows 
the branch of steady 2D axisymmetric solutions thus obtained. As the branch is followed towards 
larger L^-norm, the weakly nonlinear solutions homoclinic to A~ do indeed grow into steady circular 
fronts between A^ as postulated. 

A circular front has a well-defined radius, r = p, which may be measured by 


V{r = p) = -{V+ + V-}- 


( 2 . 1 ) 


As shown in Eig. 2(a), in the upper part the branch approaches the limit 7 = as the norm 


or equivalently the front radius p increases. This is because at p = 00, the curvature term r ^dr 
in Eq. (jl.3p vanishes and the steady circular front becomes identical to the steady ID front, which 
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Figure 2: (a) The branch of 2D axisymmetric steady states followed from the lower saddle-node at u = 5. (b) A 
sample solution profile V(x,y) at 7 = 1.844. 


exists only at 7 = 7^*^. In general, the instantaneous speed c of a circular front is a smooth function 
of 7 and p. Moreover, a circular front with large p is well described by the governing ODE in r 
restricted to an interval of width dr around r = p where 1 ^ dr <C p. Hence for 7 near and p 
large, the front speed c should depend linearly on dy = 7 — 7^*^ and the front curvature k = p~^ as 


c = c^dy + c^tv. 


( 2 . 2 ) 


The function —c{dj) ai k = 0 is plotted in Fig. 3 in Ref. [^, while d^{K) at c = 0 can be determined 
from Fig. 2(a) of this paper using Eq. (jl.9p . From the linearization of these two functions around 
the origin, we obtain numerically 


= 11.86, = -1.113. (2.3) 

In Appendix [Appendix A| we derive the analog of Eq. ()2.2p analytically for variational PDFs. 
Using the definition c = dp/dt, we can rewrite Eq. (|2.2p as the following ODE for hc: 

dn/dt = —/^^(c^dy + c,^n). (2.4) 

The nontrivial equilibrium and its eigenvalue within Eq. ( 1231 , denoted by 

Ko = (2.5) 


correspond, respectively, to the steady circular front Ao{r) and an eigenvalue in the Do spectrum 
of Ao(r). As stated above, Ao{r) becomes a steady ID front located at po = ^ oc as dy ^ 0. 

In this limit, in the Dq spectrum of Ao{r) the eigenvalue 0 is the only marginal eigenvalue and 
has multiplicity 1. The corresponding eigenfunction is the amplitude mode Aa{r) = A^ipr)^ which 
coincides with the Goldstone mode resulting from the invariance of Eq. m under translation in 
r as r ^ oo. As dy increases, the radius po of a stationary front decreases and the translational 
invariance is broken more and more strongly. As a result, the amplitude (Goldstone) mode acquires 
a nonzero growth rate as given by Eq. (j2.5p . As shown in Fig. [3l the growth rate a of the amplitude 
mode Aa (r) along the branch determined by continuation does indeed agree with Eq. (1^ for small 
dy. 


As shown in Fig. |2(a) , the LS branch first undergoes a series of ten folds, which we label as y^, 
i = 1, • • • ,10, counting from y^^. The distance between successive folds increases non-monotonically 
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Figure 3: Th e growth rate cr of the amplitude mode Aa{r) as a function of 7 along the branch of steady circular fronts 
in Fig. |2(a)] computed using numerical continuation. 



as the oscillatory tail around ^4+ is annihilated at r = 0. This behavior is reminiscent of ID collapsed 
snaking (CS), and will therefore be referred to as radially collapsed snaking (RCS). The oscillatory 
tail around ^4+ results from the fact that the spatial eigenvalues of ^4+ form a complex quartet 
denoted by (A^ ± iXi, —Xr ± iXi) with > 0. For CS, it was shown in Ref. [HI (Lemma 2.3) that 
as the front radius p increases, 7 approaches the limiting value 7 ^*^ according to 

7 — 7 ^*^ ^ exp {—2Xrp) sin {2Xip + 0 ), ( 2 . 6 ) 

where 0 is a constant. For d-dimensional axisymmetric fronts in variational PDFs, it is shown in 
Appendix [Appendix A| that if either A^^lord— 1 ^ 1 , then RCS results from a competition 
between the 0(exp(—2A^/))) tail oscillation consistent with Eq. (|2.6p and the monotonic contribution 
0{{d — l)p~^) proportional to the front curvature. If this also holds for the 1:1 FCGLE, then as p 
decreases from oc, folds should begin to appear at p = p^^^ ^ 1 , where p^^^ is given by Eq. (jA.16|) 
of Appendix [Appendix A[ In the A^ ^ 1 case, this equation can be solved asymptotically to yield 

= A;:y-logAr + log(-logAr) + 0(1)). (2.7) 

In the present parameter regime, for 7 not much larger than 7 ^*^, the spatial eigenvalues are 
A^ ^ 0.19 and A^ ^ 1.4. Hence the predicted critical radius using Eq. (|2.7p is p^^^ ^8.8, while the 
wavelength of the oscillatory tail is 27r/Ai ^ 4.5. In comparison, the measured radii at 710 and 73 
on the branch are respectively pio ~ 10.2 and pg ~ 7.5. As shown in Eig. [H the growth rate cr of 
the amplitude mode Aa{r) changes sign at 7 ^, i = 7, • • • , 10 . As a result, stable oscillons with radii 
approximately pio and pg exist, respectively, in the narrow intervals ( 79 , 710 ) and ( 77 , 7 g), and no 
stable oscillon exists with radius greater than pio. The existence of a maximum stable radius is a 
unique feature of axisymmetric oscillons in two or higher dimensions. 

As the branch reaches 70 from above, the solution becomes well localized around r = 0 with the 
front radius p comparable to the tail wavelength. Thus the RCS mechanism based on the annihilation 
of the oscillatory tail cannot account for the lower folds. The complete spectral stability of the entire 
branch is computed directly in the next section. 


2.2. Spectral stability 

Figure [4] shows the spectral stability of steady axisymmetric oscillons on the RCS branch in 
Fig. |2(a) , where 5R(cr) is the real part of the eigenvalue cr, and i is the arclength along the branch 
relative to the bifurcation point at 7 ^^. The spectra are plotted only for the lowest azimuthal 
wavenumbers m = 0 to 3 since oscillons with small radii are most likely unstable to B^ perturbations 
with small m. For each m, only the two largest values of 5R(cr), or equivalently the two most unstable 
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Figure 4: Spectral stability of steady axisymmetric oscillons on the RCS branch in Fig. |2(a)] shown in terms of the 
real part of the eigenvalue 5R((7) as a function of the arclength I measured from . For each azimuthal wavenumber 
m = 0,1,2,3, only the two most unstable eigenvalues are plotted. The zero eigenvalue for m = 1 resulting from 
translation invariance in space is not shown. The locations of the folds on the RCS branch are indicated using dashed 
vertical lines. 


or leading eigenvalues, are plotted. Thus at any the presence of two values of J?(cr) implies that 
the leading eigenvalue is real, whereas a single value of 5R(cr) implies that the leading eigenvalues 
form a complex conjugate pair. Owing to translation invariance there is always a zero eigenvalue in 
the Di spectrum, which is filtered out before the spectrum is plotted. The locations of the folds on 
the RCS branch are labeled using dashed vertical lines. 

The Do spectrum is shown in black in Fig.lH There is an amplitude mode Aa{r) whose eigenvalue 
is real and changes sign at every fold. As a result, oscillons are Dq unstable on the segments 
^SN ^ stable on the segments 71 ^ 72 and 73 ^ 74. 

Interestingly, on the segment 71 ^72, the real eigenvalue of the amplitude mode Aa{r) and a 
smaller real eigenvalue first coalesce into a complex conjugate pair, and this pair then splits back 
into the real eigenvalue of the amplitude mode Aa{r) and a smaller real eigenvalue. Thus, for 
oscillons on this segment, Dq perturbations may decay in either monotonic or oscillatory fashion 
depending on 7. This process of coalescence and splitting happens again on the segment 75 ^76, 
but this time the complex conjugate pair resulting from the coalescence of the two real eigenvalues 
becomes unstable before splitting back to two real eigenvalues, this time with the smaller eigenvalue 
corresponding to the amplitude mode Aa{r). Thus, for oscillons on this segment, Dq perturbations 
may grow or decay depending on 7 and may do so in either monotonic or oscillatory fashion. Beyond 
76, oscillons are Dq unstable on the segments 76 ^ 77, 78 ^ 79, and 710 ^ 7^*^, and Dq stable on 
the segments 77 ^ 73 and 79 ^ 710, as explained by the RCS mechanism in ^2.11 

The D^ spectra for m = 1 , 2,3 are shown, respectively, in green, red, and magenta. Since the 
odd segments j2i 722+1, '^ = 0,1, * * *, with 70 = 7^^ are always Dq unstable, we only consider D^ 
stability, m > 1, on the even segments 722+1 ^ 722+2- We observe that the leading Di eigenvalues 
are always stable, and are generally more stable than the leading D2 eigenvalues. The leading D2 
and D3 eigenvalues are stable on all even segments except 75 ^76; on part of this segment, these 
eigenvalues are real and unstable. 

In summary, oscillons on the even segments with relatively small or large radii, specifically those 
on 722+1 ^ 722+2, = 0 , 1 , 3 , 4 , are spectrally stable. However, oscillons with intermediate radii, 

specifically those on 75 ^76, are always D2 unstable and may in addition be Dq and D3 unstable. 

These instabilities may be observed by time-evolving steady oscillons with superposed small 
amplitude random noise. Figure [ 5 ] shows the time evolution at 7 = 1.857 on the segment 75 ^76, 
where the only unstable eigenvalue is 0.050 (D2). In this case, the D2 instability is triggered by the 
D2 noise in the initial condition, and the solution subsequently expands D2 symmetrically. After 
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Figure 5: Time evolution of a steady oscillon at 7 = 1.857 on the segment 75 —)■ 75 , with a D 2 noise in the initial 
condition. In this case the only unstable eigenvalue is 0.050 ( 02 )- The domain size is [—20,20] X [—20,20], and the 
snapshots are taken at t = 0,100, 200, 400. 
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Figure 6: Time evolution of a steady oscillon at 7 = 1.86 on the segment 75 —)' 76 , with a Dq noise in the initial 
condition. In this case the unstable eigenvalues are 0.037 ± O.lGz (Oo), 0.092 (D2), and 0.036 (Os). The domain size 
is [—20, 20] X [—20, 20], and the snapshots are taken at t = 0, 320, 340, 480, 540, 720. 


colliding with its images relative to the periodic domain boundaries, the solution first evolves into 
a spot homoclinic to the equilibrium ^4+; this spot then persists for a long time before eventually 
collapsing into ^4+ itself. Figure [6] shows the time evolution at 7 = 1.86 on the segment 75 ^76, 
where the unstable eigenvalues are 0.037 ± 0.16i (Bq), 0.092 (B2), and 0.036 (B3). In this case, the 
noise triggers a Bq instability, and the solution subsequently oscillates Bq symmetrically. As the 
oscillation amplitude increases, the B2 instability takes over and the subsequent evolution beginning 
with a B2 symmetric expansion is similar to Fig. [5l 

2.3. Traveling circular fronts 

As shown in ^2.11 the radius p of a steady circular front is determined locally uniquely by the 
value of the parameter 7. Traveling circular fronts are therefore expected when p is changed by dp or 
7 is changed by dy. When dp < 0 (dp > 0) and dy = 0, namely, when the front radius is contracted 
(expanded) at fixed 7, DNS shows that the circular front continues to contract (expand). Figure [3 
shows that in the latter case the circular front eventually collides with four image fronts, and evolves 
into reciprocal circular fronts centered at the corners, with A~ now embedded in A+. The solution 
eventually becomes spatially uniform when the A~ bubbles in the corners contract to zero. 

The contraction and expansion processes maintain axisymmetry provided the LS does not reach 
the domain boundary. In general, axisymmetric time evolution is described by Eq. (EH with dee = 0? 
namely 

Ai = (I + ioi) H —A + /(A). (2-8) 

Figure [8] shows radial space-time plots illustrating these processes computed from Eq. (j2.8[) . with the 
front evolution predicted using Eq. (j2.4j) shown as the dashed line. The agreement is excellent. In 
both contraction and expansion, we observe that the front speed increases as the front radius evolves 
away from its equilibrium value. The resulting front evolution is well approximated by Eq. (1221) as 
long as the front radius p remains large, or equivalently the front curvature n = p~^ remains small. 

At any given time t during the time evolution, the front radius p can be measured by Eq. m 
and the instantaneous front speed c can be computed as c = dp/dt. Moreover, DNS results in Eig.[8] 
suggest that to a good approximation, c is the instantaneous speed at which the entire front moves 
in the radial direction. Thus traveling circular fronts are well approximated by solutions to the 
radial PDE Eq. (j2.8p of the form A{r — ct) near t = 0. This ansatz turns Eq. (j2.8p into the radial 
traveling wave ODE 

— cArp = (I T iodj fdrr T ~drj A T /(A). (^*9) 
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(a) (b) 


Figure 7: Snapshots of V{x,y) for z/ = 5, 7 = 1.844, d'y = 0, dp = 1 at (a) t = 9300 and (b) t = 9700. The initial 
condition is the steady circular front in Fig. |2(b)] 




(a) 


(b) 


Figure 8 : Radial space-time plots of V(r,t) for the same parameter values as in Fig. [7| except that (a) dp = —1, (b) 
dp = 1. The dashed lines show solutions computed from Eq. 12311 for comparison. 
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Figure 9: The front speed c as a function of the front radius p at 7 = 1.844. Dotted: DNS. Solid: numerical 
continuation using Eq. iirm . The curves superpose with high accuracy. 


Using this ODE with NBC given by Eq. the steady circular front in Eig. |2(b)| with c = 0 

can be numerically continued to obtain a family of traveling circular fronts with c 7^ 0. Eigure [9] 
shows the front speed c as a function of the front radius p measured from DNS and computed using 
continuation. The agreement is excellent for the range of p over which c is measured. In agreement 
with the RCS argument leading to Eq. (EI3), c starts to oscillate as a function of p near p^^^ ^8.8. 

Eurther continuation of this branch towards more negative front speed c yields a family of 
contracting pulse rings. A contracting pulse ring is essentially a ID stationary pulse bent into 
a circle in 2D, and thus its speed c is proportional to its curvature /^; see Appendix [Appendix A| 
for an analytical description of contracting/expanding pulse rings in variational PDEs. As shown 
in Eig. [TOlab the proportionality constant = c//^ is negative. The governing ODE for the ring 
radius, ^ = c = can be solved to yield 

= ( 2 . 10 ) 

with the constant to determined by the initial condition. As shown in Eig. flOl b). the contracting 
pulse ring taken from this branch is amplitude-unstable and thus rapidly evolves into a thicker 
contracting pulse ring. The latter contracts to zero in a parabolic fashion as predicted by Eq. (I2TQ1) 
with < 0; nontrivial dynamics are visible just prior to extinction when the radius of the ring 
becomes comparable to the pulse width. We remark that the evolution of a ID stationary pulse 
bent into a closed 2D curve of an arbitrary shape may be analyzed following previous work on the 
Swift-Hohenberg equation [l^ . 


3. Localized target patterns 

The properties of time-independent Type-II LS of Eq. (HU in ID are explored in Refs. Ulii. 
These structures consist of a Turing pattern bifurcating from A+ embedded in the homogeneous 
state A~. The associated bifurcation diagram consists of a single solution branch oscillating 
between two distinct limiting parameter values; the oscillations reflect the growth of the Turing 
pattern via roll insertion at the center This bifurcation behavior is termed defect-mediated 

snaking (DMS), the limiting parameter values are called the snaking limits, and the parameter 
interval inbetween is called the snaking interval. At fixed u = 7 snaking takes place in the interval 
^^DMS^^DMS] _ [2.8949,2.8970]. Outside the snaking interval, Type-II LS depin by creating or 
destroying rolls in the interior of the wavetrain via successive phase slips [2^. These phase slips 
occur at the center of the LS when 7 is near either limit of the snaking region, but off-center for 7 
farther away. 
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Figure 10: A contracting pulse ring at 7 = 1.844. (a) The speed c as a function of the radius p. (b) Radial space-time 
plot of R(r, t) with initial condition corresponding to c = —0.04. 


In this section we study a family of 2D localized target solutions of Eq. (jl.ip formed by embedding 
into A~ an axisymmetric Turing pattern that bifurcates supercritically from ^4+ (Fig. [I{a)). As 
in Refs. [13, [HI we fix z/ = 7 and vary 7 as the bifurcation parameter. In ^3.11 we compute the 
bifurcation diagram of steady localized target patterns, and compare their bifurcation behavior with 
ID DMS and 2D RCS. In ^3.2l we compute the spectral stability of solutions along this branch, and 
identify instabilities with possible relevance to experiments. In ^3.31 we compute the radial depinning 
of localized target patterns, and compare their temporal dynamics with the depinning of Type-II 
LS in ID. 


S.l. Core-mediated snaking 

The branch of 2D axisymmetric steady states similar to Fig. 2(a)| but computed at z/ = 7 is 
shown in Fig. ll(a)[ Overall this branch is similar to the 2D RCS branch shown in Fig. |2(a)[ 
where the curvature of the circular front causes the upper part to converge monotonically to a 
single limiting parameter value. In contrast to Fig. |2(a)[ in the present case the branch undergoes 
^-shaped snaking on narrow intervals of 7. The separation between successive intervals decreases 
as the L^-norm increases and the solutions high up the branch take the form of localized target 
patterns (Fig. |ll(b)[ ). As shown in the video accompanying Fig. [TTl the localized target patterns 
in Eq. grow through a central defect at the core r = 0, similar to ID DMS. This bifurcation 
structure is hereafter referred to as core-mediated snaking (CMS). 

The snaking mechanism of localized target patterns in the 1:1 FCGLE may be compared with the 
behavior of the corresponding structures in SH23 [HI . When d = 1 new rolls in SH23 are nucleated 
at the location of the fronts while in the 1:1 FCGLE they are inserted by the division of the central 
“cell” or “defect”. When d > 1 localized target patterns in both equations possess a spike at r = 0. 
As a result, as d increases, the snaking mechanism in SH23 changes from ID standard homoclinic 
snaking to 2D CMS via a stack of isolas. In contrast, in the 1:1 FCGLE the defect and the spike are 
both located at r = 0 and as d increases the snaking mechanism is expected to transition continuously 
from ID DMS to 2D ^-shaped GMS without the break-up of the snaking curve. We remark that in 
either SH23 or the 1:1 FGGLE, since the core plays the role of a “pacemaker” that emits new rings, 
the solutions with a local maximum/minimum at r = 0 (known as type A/B solutions in Ref. |26| ) 
necessarily exist on the CMS branch in an alternative fashion. 

To test this conjecture, we have continued solutions to Eq. dm in the dimension d; the result 
of this continuation at 7 = 2.897 is shown in Fig. [121 In order that the snaking in d corresponds to 
snaking in 7 at nearby values of d, we plot the solution branch in the (—d, ||A||a,) plane, such that 
the sharper folds are located on the right. In this way a pair of adjacent C-shaped segments forms 
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Figure 11: (a) The branch of 2D axisymmetric steady states ai u = 7. A video showing the growth of the LS as the 
branch is followed from the lower saddle-node can be viewed at this URL (b) A sample solution profile V{x^y). 


the body of the letter Near d = 1 the snaking in d takes place mostly in d > 1 and resembles 
DMS except that the folds corresponding to maxima and minima in V{r) at r = 0 no longer align 
with each other. As d increases a hysteresis loop appears on the bottom segment of every pair of 
adjacent C-shaped segments, while the snaking interval in d becomes narrower. At d = 2, the two 
additional segments forming this hysteresis loop can be regarded as the (short) tail of the current 
and the (long) head of the preceding letter Thus one extra ring is added at the core every time the 
branch passes through a ^-shaped segment. At large d (e.g. d = 3), the body of the ^ disappears via 
two hysteresis bifurcations, although the head and tail remain. We have not continued the branch 
to even larger d and thus cannot ascertain whether another hysteresis bifurcation would eventually 
lead to a bifurcation curve without folds, but from Fig.[T2]the branch appears to approach a straight 
line in the (—d, ||A||a,) plane as d ^ oo. 

Apart from continuation in d for 7 = 2.897 near we have also examined other slices at 

constant 7 but these continuations do not lead to axisymmetric steady states in higher dimensions. 
Specifically, for smaller 7 (near 7^^^*^) the branch snakes towards d < 1 while for intermediate 7 
(about halfway between and branch snakes around d = 1. 


3.2. Spectral stability 

Figure [13] shows the spectral stability of steady axisymmetric oscillons on the CMS branch in 
Fig. ll(a)[ plotted using the same conventions as in Fig. jH As in Fig. IH oscillons on even segments 
with relatively small radii are spectrally stable. In contrast to Fig. IH oscillons with intermediate 
radii may be unstable not only for m = 0, 2, 3, but also for m = 1. Also, oscillons with relatively 
large radii, which resemble localized target patterns, are generally B^ unstable for m > 1 due to 
azimuthal instability of the rings. 

Figure (in shows the time evolution of a steady oscillon at 7 = 2.922 on the segment 79 ^ 710, 
where the unstable eigenvalues are 0.022 ±0.12i (Bq) and 0.076 (Bi). In this case, the noise initially 
excites a Bi instability, and the solution subsequently develops a Bi symmetric interior. The Bq 
symmetry is subsequently restored, triggering the Bq instability. As a result, the solution first 
oscillates and then expands Bq symmetrically, before colliding with its images and eventually evolving 
into a hexagonal pattern. 

For m > 1, the locations on the RCS or CMS branch where the B^ stability changes correspond 
to bifurcation points to steady m-fold symmetric oscillons. Further from the onset, these m-fold 
symmetric oscillons may become stable, in which case they may be identified with the m-fold 
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Figure 12: The branch of d-dimensional axisymmetric steady states followed in d from a ID localized state (d = 1) 
on the DMS branch at i/ = 7, 7 = 2.897. 



Figure 13: Spectral stability of steady axisymmetric oscillons on the CMS branch in Fig. |ll(a)[ plotted using the 
same convention as in Fig. [4| For clarity only the first ten folds are indicated using dashed vertical lines. 



Figure 14: Time evolution of a steady oscillon at 7 = 2.922 on the segment 79 —)► 710 , with a Di noise in the 
initial condition. In this case the unstable eigenvalues are 0.022 ± 0.122 (Dq) and 0.076 (Oi). The domain size is 
[—20, 20] X [—20, 20], and the snapshots are taken at t = 0, 70, 90,106,110, 800. 
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symmetric oscillon “molecules” observed in some experiments HSUS. Thus, dipoles formed by 
two out-of-phase oscillons may bifurcate from a Bi bifurcation point, dimers formed by two in-phase 
oscillons may bifurcate from a D2 bifurcation point, and so on. A detailed study of the interaction 
between simple oscillons and the bifurcation of oscillon molecules is left for future work. 


3.3. Radial depinning and traveling circular fronts 

Owing to the presence of azimuthal instabilities, fully 2D DNS of localized target patterns 
generally produce fully 2D localized states, as discussed further in ^4.11 On the other hand, azimuthal 
instabilities may be artificially suppressed by imposing axisymmetry on the solutions. This procedure 
allows us to time evolve localized target patterns using Eq. (12.Sp and to study radial depinning of 
the bounding front, as done in ^2.31 Since the front speed c and the dynamically selected interior 
wavenumber k are in general coupled [IBl we assume in the following that the wavenumber selection 
is largely passive (i.e., that c depends only weakly on k). In this case the speed c may be regarded as 
a function of the control parameter 7 and the front radius p only, as in ^2.31 Hence to obtain c 7^ 0, 


we can start with a localized target pattern at a particular fold on the CMS branch in Fig. |ll(a) 
and either change the radius p by dp or the parameter 7 by dy; the outcome is similar. Figure [15 
shows radial space-time plots computed for different choices of the fold location and perturbation 
dy. Figures [15 (a) and |15(b')] show the time evolution near and far from the CMS region, respectively, 
in both cases starting from an initial condition at a left fold on the CMS branch; Figs. 15(c) and 
|15(d)| show the corresponding results for an initial condition at a right fold instead. In all four cases, 
the front speed increases as the front radius evolves away from its equilibrium value, and individual 
rings are destroyed or created in the interior of the structure via repeated phase slips. Thus, radial 
depinning of localized target patterns inherits the key features of both traveling circular fronts and 
depinning of Type-II LS in ID. A particularly interesting new feature is that near the CMS region 
the phase slips occur right at the core, while far from the CMS region the phase slips occur both 
right at the core and at a nonzero radius away from the core. 

Of particular interest is the behavior of the core at r = 0. To the left of the left fold the structure 
contracts, but the way it does so depends on whether the perturbation is small in which case the 
frequency with which phase slips occur in the core accelerates as the target contracts (Fig. 15(a)[ ), 
or the perturbation is large in which case the majority of the phase slips take place away from the 
core (Fig. |15(b)[ ). Similar behavior is seen when the target expands. If the perturbation is small 
the phase slip frequency in the core continues to increase as the target expands (Fig. |15(c)[ ), but it 
decreases with time when the perturbation is large (Fig. |15(d)[ ). Thus the different magnitudes of 
the perturbation have distinct consequences. A small perturbation destabilizes the core ring first 
since curvature narrows the range of stable rings. The resulting phase slip eliminates a ring at 
the core but the inward front motion triggers additional phase slips at the core. As the structure 
contracts more and more rapidly the phase slip frequency necessarily increases. In contrast, when 
the structure contracts because of a large perturbation the inward motion of the front compresses 
mostly the outer rings (the curvature provides the inner rings with increased stiffness) leading to 
phase slips in the vicinity of the front. The converse happens when the structure expands. Thus the 
system “remembers” the initial perturbation that led to the contraction or expansion, in contrast 
to the ID case in which the phase slips near the fronts on either side gradually migrate inwards 
to or outward from the center x = 0 as the structure contracts or expands. We attribute these 
differences to the reduced role played by curvature in the outer regions of the pattern and hence to 
radial inhomogeneity of the contraction or expansion processes. 

As in the ID case [HI, Eq. (|2.9[) contains codimension-1 traveling circular fronts between 
for either sign of the front speed c. Figure n a) shows two branches of contracting/expanding 
circular fronts computed at c = —0.05 (blue) and c = 0.05 (red), which are found to exist on 
opposite sides of the ^-shaped CMS branch of steady localized target patterns (black). Near c = 0, 
these two branches are expected to approach the left and right snaking limits of the CMS branch 
and thus remain separated by a finite distance. The origin of this phenomenon is the so-called 
“reversibility-breaking” bifurcation on A+ at c = 0 [25|. As c crosses zero, the two pairs of spatial 
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Figure 15: Radial space-time plots of V{r,t) at = 7 and near the fold at (a,b) 7 = 2.897441, (c,d) 7 = 2.897556, 
with the perturbations (a) ^7 = — 1 x 10 “"^; (b) ^7 = — 1 x 10 “^; (c) d'j = 5 x 10 “^; (d) d'j = 2 x 10 “^. 
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Figure 16: (a) Two branches of traveling circular fronts computed at c = —0.05 (blue) and c = 0.05 (red), and 
the ^-shaped CMS branch of steady localized targets (black), (b) The bifurcation parameter 7 as a function of the 
curvature k ai c = —0.05. (c) The logarithm of the core amplitude \A(r = 0) — A+| as a function of the front radius 
p at c = —0.05. 


eigenvalues and =bi/c2 {ki ^2 > 0) of ^4+ cross the imaginary axis simultaneously in opposite 

directions. As a result, the unstable eigendirections of A+ change by a finite amount, and so does 
the value of 7 at which the heteroclinic orbit between exists. 

For 0 < |c| ^ 1, the real part of the unstable eigenvalue of A+ scales as 0(c) [IBl. Thus as 
discussed in 92.11 the branch of traveling circular fronts between should exhibit RCS. FigurefTbfb) 
shows the bifurcation parameter 7 as a function of the curvature k at c = —0.05. The solid dot 
labels the critical radius = 51 beyond which folds disappear. At this point one calculates that 

Xr = 0.087, so the predicted critical radius based on Eq. (loi is = 39. Thus even for relatively 
small Xr^ Eq. (1^ underestimates the critical radius p ^^^. To understand this discrepancy, we note 
that the underlying assumption of Eq. (Ell) that the tail oscillation is 0(exp(—2Ar-p)) presumes that 
the core amplitude is 0(exp(—A^p)). While this holds true for ID LS, 2D axisymmetric LS are 


often found to have larger amplitude cores [21 


Eigure He) shows the logarithm of the 


amplitude \A{r = 0) — A+| as a function of the front radius p at c = —0.05. It can be seen that the 
decaying exponent is indeed A^* for large p, but becomes less than A^ for smaller p including p values 
around p^^^. 
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4. Localized hexagonal patterns 

In this section we study a family of 2D localized hexagonal patterns consisting of fully 2D 
hexagonal patterns embedded in A~. This family is further divided into circular ( ^4.ip and planar 
(ESI localized hexagonal patterns. We distinguish between two cases: the supercritical case 7 < 7^ 
and the subcritical case 7^ < 7 < 7^. As in 0 we fix = 7 and vary 7; the Turing bifurcation 
point is 7^ = 3.0845. 


4 . 1 . Circular localized hexagonal patterns 

In the supercritical case fully 2D LS can be found by time-evolving localized target patterns 
in 2D. The rings rapidly break up into hexagons as a result of azimuthal instabilities triggered by 
the discretization. Figure [T3 shows snapshots of the initial stages of this 2D time evolution at 
7 = 2.8989 with the localized target pattern in Fig. |ll(b)| used as initial condition. The fully 2D 
pattern that emerges consists of an approximately hexagonal pattern filling the interior of a circular 
front. Since the 2D solution is discretized along the x and y directions, the initial pattern contains 
multiple defects arranged into a B4 symmetric configuration (Fig. 17(a)[ ). Shortly thereafter, this 
B4 symmetric pattern breaks up from the center and the defects reorganize into a pair of parabolas 
with B2 symmetry (Fig. |l7(b)[ ). This symmetry breaking process is reminiscent of the creation of 
B2 symmetric dimers or B3 symmetric trimers from a Be symmetric hexagonal lattice of oscillons 
observed in experiments [30|, l20| . 

The subsequent evolution of the state in Fig.flTla) is driven by the presence of the defects. These 
take the form of the penta-hepta defects familiar from other haxagon-forming systems (these are 
located along the x-y axes in the figure and nearer the center); however, free pentagonal defects 
are found near the periphery of the structure where curvature effects are significant. The latter 
serve as loci of phase slips whereby new cells are inserted (or deleted) from the structure leading 
to readjustment of nearby cells that in turn triggers the formation of new defects. As a result the 
whole structure undergoes slow, phase-slip driven evolution. Figure fTTI b) shows a later state of the 
system: the structure has almost the same size and shape but the defect distribution differs, with the 
defects forming prominent chains at later times (see the video at this URL). In contrast to ^2.31 and 
^3.31 this defect-driven evolution cannot be axisymmetric, and so the front bounding the hexagons 
does not remain perfectly circular at all times. Simulations indicate that there is an approximate 
critical radius po such that the bounding front contracts (expands) when the front radius p < po 
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(a) 


(b) 


Figure 18: Snapshots of V{x,y) at u = 7, and (a) 7 = 2.8978, t = 800; (b) 7 = 2.8989, t = 4080, showing respectively 
a contracting and an expanding circular localized hexagon. Videos of the evolution can be viewed at this URL, and 
this URL 


{p > po). When p is far from po, the azimuthally-averaged front speed c depends mainly on 7 and 
the azimuthally-averaged front radius p evolves in a manner similar to the circular traveling fronts 
described in ^2.31 However, when p is near po, the average front speed c becomes more sensitive to 
the details of the hexagonal pattern, which makes it difficult to determine po precisely or stabilize 
the evolution by feedback control. 

When p is near po, the bounding front expands (contracts) slowly and hexagons can be created 
(annihilated) anywhere within the region they occupy. When p is further from po, the bounding front 
expands (contracts) more rapidly and hexagons tend to be created (annihilated) at a fixed distance 
from the bounding front. This faster evolution is analogous to off-center depinning of Type-II LS in 
ID. For both slow and fast evolution, the bounding front remains almost perfectly circular. However, 
when p is far from po, fast viscous shocks may propagate azimuthally along the bounding front and 
these create (annihilate) hexagons on a yet faster timescale. This evolution has no ID analog and 
is referred to as ultra-fast. Figure |18(a) shows a snapshot of a contracting solution obtained from 
an initial p that is slightly below po; in this case the initial evolution is slow but transitions to 
fast and subsequently ultra-fast contraction once p becomes sufficiently small. Figure p~8(b)| shows 
a snapshot of an expanding solution for which the initial p is slightly above po; in this case once 
fast expansion sets, the circular localized hexagonal pattern collides with its images, and evolves 
into a reciprocal circular localized hexagonal pattern centered on the corners. In contrast to the 
contraction case, ultra-fast expansion does not take place when p becomes very large. Once the 
reciprocal circular localized hexagonal pattern forms, hexagons are always created right behind the 
bounding front; of course this front has negative curvature, in contrast to the bounding front in 
Fig. |18(a)| The domain is ultimately filled with a hexagonal pattern containing multiple chains of 
penta-hepta defects inherited from the initial evolution. 

The presence of front pinning in localized target patterns and the apparent absence of pinning 
in circular localized hexagonal patterns can be understood by comparing the orientations of the 
bounding front and the periodic pattern within. In a localized target pattern, the circular front 
is always pinned to a periodic pattern in the direction perpendicular to itself, i.e., in the radial 
direction. However, in a circular localized hexagonal pattern, the circular front cannot be pinned to 
a hexagonal pattern that is only periodic along three directions. 
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^.2. Interpretation: Single and double front strueture 

As illustrated in Fig. [I{b), in 2D the Turing bifurcation on A+ results generically in stripes and 
hexagons. The bifurcation to the resulting unstable branch of stripes S is supercritical as shown in 
Ref. [ 2 ^. The bifurcation to hexagons is transcritical, with the bifurcating hexagons unstable on 
both sides of the bifurcation. As detailed next, the up-hexagons consisting of a hexagonal array 
of peaks (in terms of V = ^(A)) bifurcate subcritically, i.e., into 7 > 7^, while the down-hexagons 
H~ consisting of a hexagonal array of dips bifurcate supercritically, i.e., into 7 < 7^. The former 
turn around in a saddle-node bifurcation at 7 = 7^ and acquire stability on the upper part of the 
branch (Fig. [Ifb)). Thus in 7^ < 7 < 7^ the stable state coexists with the stable A+ 
state, while for 7 < 7^ only the state is (locally) stable and so coexists with the stable A~ 
state instead. 

In the subcritical regime 7^ < 7 < 7^ stationary fronts exist between and A+, and these 
are a consequence of the pinning of the front to the hexagonal structure on one side. As illustrated 
in Fig. I^a.b). these stationary fronts can be assembled into steady fully and planar localized 
hexagonal patterns. Figure [ 19 ] shows examples of both types computed via DNS with the initial 
condition taken to be a hexagonal (rectangular) patch of hexagons embedded in A+. The 
hexagons are obtained by incrementing 7 adiabatically from 7 < 7^ past 7^ and into 7 > 7^. 
The computations are done with the horizontal and vertical domain sizes Lx = 200\/3 and Ly = 200 
for which the wavenumber k is close to the critical Turing wavenumber at 7 ^, and Ly is occupied 
by exactly A" = 49 rows of hexagons. The resulting hexagons remain amplitude stable for 7 < 7^, 
where 7^ ^ 3.18 and are well approximated by 

A = A~^ + A(cos (kx) -h cos {k{x -h V3y)/2) -h cos {k{x — V3y)/2)) ( 4 . 1 ) 

with A+ = 1.044 — 0 . 495 i and A = — 0.004 -h 0 . 116 i at 7 = 3.14 > 7^. These hexagons are thus 
up-hexagons (i.e., V = ^{A) > 0 ); indeed as seen in Figs. [T 3 & [191 ^ consists of 

a hexagonal array of peaks. We remark that for some parameters the wavenumber k is found to 
be Eckhaus unstable. For both fully and planar localized hexagonal patterns, the front between 
and A+ selects a preferred hexagonal wavelength, and the subsequent dynamics depend on 
the choice of 7. In the supercritical regime 7 < 7^ the stable invades the unstable A+ via a 
pulled front. As 7 increases past 7^, the pulled front transitions into a pushed front between the 
bistable and A+. In the subcritical regime 7 > 7^, the pushed front expands for smaller 7 and 
contracts for larger 7; videos of the expansion and contraction of fully localized hexagonal patterns 
(Fig. EOl^a)) at 7 = 3.1438 and 7 = 3.1441 can be viewed at this URL and this URL, while videos 
of the expansion and contraction of the related planar localized hexagonal patterns (Fig. EOj^b)) at 
7 = 3.1555 and 7 = 3.1557 can be viewed at this URL and this URL, Interestingly, in the later 
stages of expansion the hexagonal pattern has line defects along the boundaries and thus does not 
fill the domain perfectly. For intermediate 7 the front is expected to be pinned to the hexagonal 
pattern over a finite interval of 7 known as the pinning region. The above DNS results suggest that 
the pinning region may be very narrow, which is consistent with the slowly varying front profile and 
the proximity to the Turing bifurcation 0 ; Fig. [TCTa) and Fig. [T 9 lb) show almost steady fully and 
planar localized hexagonal patterns at 7 = 3.14405 and 7 = 3 . 1556 . Similar solutions have been 
studied in detail using numerical continuation for SH 23 [22|, while associated depinning dynamics 
present outside the pinning region have been studied numerically in the Lengyel-Epstein model [l 5 | 
and in a crystallization model by Archer et al. Q. 

Theory 0 shows that in the subcritical regime 7^ < 7 < 7^ the front between and A+ 
is a single pushed front. However, in the supercritical regime 7 < 7^ the front between and 
A+ generally exhibits a double front structure, with a first front connecting the stable state 
to the unstable S state, and a second front connecting the S state to the unstable A+ state. In 
Ref. [H the transition between the single and double fronts near 7^ has been studied in some detail 
using appropriate amplitude equations 00. One finds that sufficiently far from 7^ the speed of 
the S-A'^ front is determined by the marginal stability criterion as appropriate for a pulled front 
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(b) 


Figure 19: Snapshots of V{x,y) at (a) 7 = 3.14405, t = 5000; (b) 7 = 3.1556, t = 10000, showing almost steady fully 
and planar localized hexagonal patterns in the subcritical regime. 


0 and likewise for the H^-S front. However, as the distance 7^ — 7 is reduced, the H^-S speed 
departs from the predictions of the marginal stability theory and is determined instead nonlinearly. 
This is a consequence of an orbit-flip whereby the front approaches S with a faster decay rate than 
expected generically. This transition leads to a speed that connects smoothly to the pushed front 
speed computed in the subcritical regime. 

The circular localized hexagonal patterns described in are supercritical localized hexagons 
with embedded in A~ (see Fig. nil). Since in this parameter regime is unstable and coexists 
with stable A ~, the circular fronts that form between and A~ are strongly nonlinear and cannot 
be described by weakly nonlinear theory in Refs. [30. We conjecture, however, following [H, that 
in this case the front between and A~ exhibits an analogous double front structure, with a first 
front connecting the stable state to the unstable state S', and a second front connecting the S 
state to the stable A~ state; see Fig. EIc) for an illustration. The stripe state S is visible in the 
transition region between these two fronts as a ring enclosing the hexagonal state. The S-A~ front 
is described by Eq. (j2.8|) and studied in ^3.31 its speed may be inward or outward depending on 7 
and the front radius p. Meanwhile, the H^-S speed is always outward since the stable state 
necessarily invades the unstable S state, and always larger than the S-A~ speed. As a result, a 
single stationary stripe is observed along the boundary of the circular localized hexagonal patterns 
and its speed is expected to be slaved to the slower contraction or expansion speed of the hexagonal 
state driven by phase slips within. 

The above notions become clearer in the case of localized hexagonal patterns organized in a 
stripe-like structure, referred to as planar localized hexagonal patterns^ as illustrated in Fig. [2Qr d) 
and discussed next. 


4.3. Planar localized hexagonal patterns 

Planar 2D LS can be found by time-evolving localized stripe patterns (Fig. 21(a)[ ) in 2D, with 
superposed small-amplitude noise to excite transverse instabilities of the stripes. Figure (b) | shows 
snapshots of the initial stages of this 2D time evolution. The fully 2D pattern that emerges consists 
of hexagons that align themselves along the stripes with possibly multiple defects. The subsequent 
dynamics of this hexagonal pattern in the A~ background depend on the value of 7 used in the time 
evolution. However, in each case the fronts connecting the structure to A~ remain approximately 
straight and hence the double-front picture sketched above is no longer complicated by curvature 
effects. 
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Figure 20: Sketch of subcritical and supercritical localized hexagons, which respectively exist for 7 ^ < 7 < 7 ^ 
and 7 < 7 ^ in Fig. (20). (a) Subcritical fully localized hexagons, (b) Subcritical planar localized hexagons, (c) 
Supercritical circular localized hexagons, (d) Supercritical planar localized hexagons. 



(a) (b) 


Figure 21 : Snapshot of V{x,y) at (b) t = 80, starting from localized stripes at (a) t = 0. Parameters: i/ = 7, 
7 = 2.8972. 


22 
















































(a) 


(b) 


Figure 22 : Snapshots of V{x,y) at u = 7, and (a) 7 = 2.8955, t = 400, (b) 7 = 2.899, t = 1500 showing, respectively, 
the contraction and expansion of planar localized hexagons. Videos of the evolution can be viewed at this URL and 
this URL The shock-like structures propagating rapidly along the fronts are easily identified. 


For smaller 7, the planar localized hexagonal pattern eventually contracts into a planar symmetric 
pulse, namely a ID symmetric pulse invariant in the y direction. In this process the hexagonal cells 
are repeatedly annihilated by viscous shocks traveling along the bounding fronts, similar to the 
ultra-fast contraction phase of circular localized hexagonal patterns. A snapshot of the resulting 
state is shown in Fig. | 22 (a) at 7 = 2 . 8955 . For larger 7, the planar localized hexagonal pattern 
expands into a possibly imperfect hexagonal pattern that fills the 2 D domain. In this process new 
hexagonal cells are created in the interior by phase slips centered on two lines of defects that result 
from the slow separation of the bounding fronts, similar to the fast expansion phase of circular 
localized hexagonal patterns. A snapshot of the resulting state is shown in Fig. | 22 (b)| at 7 = 2 . 899 . 

Despite the outward (inward) motion of the bounding fronts at larger (smaller) 7, we have not 
found steady planar localized hexagonal patterns at any intermediate 7. As shown in Fig. 23 (a)[ 
at 7 = 2.8972 the boundary shocks and interior defects compete with each other, leading to 
non-unidirectional front motion. As a result the final state (Fig. | 23 (b)[ ) is considerably displaced 
from the initial centerline {x = 0) owing to the stochastic creation and annihilation of hexagonal 
cells. More remarkably, this final state provides an example of a quasi-steady localized state with 
a complicated inner structure, which is quite unexpected in spatially homogeneous pattern forming 
PDFs. The inner structure consists of a few rows of hexagonal cells interspersed with defects arranged 
in a frustrated configuration that overall neither contracts nor expands, although the competition 
between these two tendencies manifests itself occasionally as local pulsations. A remarkable video 
(see this URL) shows the evolution of this state. The earlier parts of the video show clearly the 
propagation of phase slips, whereby a phase slip triggers an adjacent phase slip etc. The later parts 
of this very long simulation show persistent dynamics in an apparently equilibrium state, i.e., a state 
of a more-or-less fixed width. This residual time dependence may be associated with the fact that 
the state is invading the unstable stripe state S from the inside while the stable background state 
A~ invades from the outside, just as in the case of the almost steady circular localized hexagonal 
patterns. The resulting fronts are approximately stationary but select wavelengths dynamically; 
these wavelengths may lie outside the stability balloon of the hexagonal state leading to phase 
slips and may, moreover, be incommensurate with the width of the structure, thereby providing an 
additional source of disequilibrium. Analogous planar localized hexagonal patterns in the subcritical 
case are shown in Fig. fT^ b). 

The difference between the dynamics of circular and planar localized hexagonal patterns may be 
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(a) (b) 

Figure 23: Snapshots of V{x,y) at (a) t = 10000 and (b) t = 30000 showing the competition between contraction and 
expansion of planar localized hexagons. Parameters: u = 7, 'y = 2.8972. A video can be viewed at this URL 


attributed to the effective boundary conditions experienced by the hexagons. In the circular case, 
the bounding front is a single closed curve that acts as a free boundary for the hexagons. As a result, 
the wavelength of the hexagons can change freely in response to front motion. However, in the planar 
case, the two bounding fronts are localized in x and periodic in y, and thus act as free boundaries 
for the hexagons only in the x-direction. Thus the x-wavelength of the hexagons can change freely 
in response to front motion, but the ^-wavelength is fixed by the periodic boundary conditions in y. 
This anisotropy of the hexagonal pattern may help explain the coexistence of contraction associated 
with the fronts and expansion from the interior of the structure. 


5. Conclusion 

In this paper we have characterized 2D LS in the 1:1 FCGLE with different inner patterns and 
front shapes. The parameters are chosen near a codimension-two point between a steady front 
between two spatially homogeneous equilibria and a supercritical Turing bifurcation on one of them. 
Steady circular fronts in the Turing-stable regime and localized target patterns in the Turing-unstable 
regime respectively exhibit new bifurcation structures termed radially collapsed snaking (RCS) and 
core-mediated snaking (CMS). Axisymmetric oscillons on the RCS and CMS branches are found 
to be spectrally stable over a wide parameter interval, and subject to various types of radial and 
azimuthal instabilities otherwise. New temporal dynamics involving front contraction or expansion 
are found using radial and fully 2D time evolution, the latter leading to fully 2D LS including 
supercritical circular and planar localized hexagonal patterns. 

The stable axisymmetric oscillons in the 1:1 FCGLE may be compared with the 2:1 FCGLE [^ . 
In the 2:1 case, stable 2D oscillons exist on a much narrower parameter interval than stable ID 
oscillons. In contrast, in the 1:1 case, stable 2D oscillons exist on multiple segments of the RCS 
and CMS branches, and the width of the existence interval is comparable to stable ID oscillons. 
The existence of multiple folds on these branches results from the existence of a supercritical Turing 
bifurcation near a stable front, which is not known to occur in the 2:1 case. However, stable 2D 
oscillons may be similarly present in other harmonically forced systems relevant to experiments. 

The subcritical localized hexagonal patterns in the 1:1 FCGLE (Fig. [T^ may be compared 
with SH23 [m. In particular, the bifurcation diagrams of steady solutions may be computable 
using numerical continuation following Ref. [22|. In the supercritical case new types of temporal 
dynamics include the slow expansion (contraction) of the hexagonal pattern driven by phase slips. 
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and the non-unidirectional front motion of planar localized hexagonal patterns near the equivalent of 
a Maxwell point. It is possible that in large domains the occurrence of phase slips may be treated as 
a stochastic process, with a corresponding stochastic description of the front motion. In particular, 
it may be that on particular timescales, the front motion near the Maxwell point can be viewed as a 
random walk. These features may be further understood by locating the penta-hepta defects through 
computing the coordination numbers of the individual hexagons, and studying their lifetimes. These 
ideas are left for future work. 

We remark, finally, that the present system also exhibits parameter regimes where the 
homogeneous equilibria and the hexagonal pattern are all simultaneously stable. In 
particular, the “Maxwell” point between may he inside the pinning region between and 
A'^. Near this tristable situation, the front speed competition may produce additional interesting 
dynamics. 
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Appendix A. Axisymmetric fronts in variational PDEs 

In this Appendix we consider the following class of variational PDEs for G M: 

Ut = C[V'^]u + £[V^] = Af{u) = ^ (A.l) 

i i 


and seek d-dimensional axisymmetric traveling wave solutions of the form u{r — ct), where r > 0 
denotes the radial coordinate. At t = 0, these solutions satisfy the following ODE 


— cUr = jC 



d — 1 
r 


dr 




(A.2) 


To derive an expression for the instantaneous speed c, we expand the operator C as 


C 



d-l 

r 


dr 


jC[dr 


-drC'[dr 


o 


multiply Eq. (IA.2P by Ur and integrate over r to obtain 


(A.3) 


d — 1 /d— 

u^dr = — J Ur {C[drr]uAf{u)) dr — J Ur - drC'[drr\udr ^ O 


Eor the integral of to be bounded, we require Ur{r oo) = 0 or equivalently u{r oc) = 
where u'^ is a constant. 

The first term on the right side of Eq. n becomes 


M = H(r = 0) - H(r = oo), (A.5) 

where the function H (r) denotes the following indefinite integral 

H{r) = J Ur {jC[drr]u ^ dr. (A-6) 

The integrand takes the form of a total derivative. After integration, the function H{r) depends 
only on the j-th derivatives u^^\r)^ j = 0,1,2, • • •, and can be identified as the ID steady state 
Hamiltonian of Eq. (EH). 
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Now let us assume that Ur nearly vanishes except in an 0(1) interval around r = p for p ^ 1 , 
or equivalently that u{r) is an axisymmetric front with large radius p connecting two flat states 
u{r 0) = u~ and u{r oo) = u'^. In this case we can simplify the second term on the right 
side in Eq. (Toil by approximating r ^ with p ^ in the integrand and then integrating by parts to 
obtain 

d-1 

1C = - / UrrC.'[drr\udr^ 

P Jo 

with an error term 0{{d — l)p~^) that can be combined with the third term of Eq. (jA.4p . 

In summary, we obtain from Eq. (IA.4I) the following analytic expression for the front speed c 

c=^y* u^dr^ {M + /C) + O ^ 

where JH and 1C are given by Eqs. and (TOt . Since for large p the r ^dr term in Eq. (jA. 2 p 

can be omitted without affecting the relevant dynamics to leading order, the integrals in Eqs. (jA.7|) 
- (|A. 8 p can be approximated by replacing u{r) with the ID front profile at the same parameters. 
Physically the front motion is driven by both the free energy (or equivalently Hamiltonian) difference 
M. and the surface (or inter facial) energy 1C. To leading order At evaluates to 

Mo = H{u-)-H{u-^), (A.9) 

and may take two qualitatively different forms: 

• The two equilibria are equivalent when u{r) is a pulse {u~ = u'^) or an Ising front 
{u~ = —14+ and M{—u) = —M{u)). In this case 

Mo = 0 . (A.IO) 

• The two equilibria i 4 + are inequivalent when u{r) is a Bloch front. In this case Mo = ^ only 
at the Maxwell point 7 = 7 m with 7 denoting the control parameter. Near 7 m 

Mo (XJ -'JM- (A.II) 

The term 1C is proportional to the front curvature k. = p~^. It follows from Eq. (jA.Sp and the 
definition c = dp/dt that steady axisymmetric fronts are unstable with respect to radius change 
when IC/K. < 0, and stable otherwise. In particular, for one-component reaction-diffusion equations, 
or equivalently when >C[V^] = V^, these fronts are always unstable because 

POO 

lC/K = —{d—l) / u^drKO. (A. 12 ) 

Jo 

However, when higher order terms like are present in >C[V^], as is the case for the Swift-Hohenberg 
equation and its generalizations, these fronts may be stable in certain parameter regimes since 1C/k 
is an alternating sum that can take either sign. 

Appendix A.l. Radially collapsed snaking in (I -|-e) dimensions 

The function c(p) determines the front evolution via dp/dt = c(p); this function in turn depends 
on M and 1C. The term M depends on p due to the deviation of u{r) from u~ near r = 0 , i.e. the local 
front profile. This profile is monotonic (oscillatory) in r if the slowest unstable eigenvalue(s) of u~ 
are real (complex). Hereafter we consider only the complex case, which occurs when >C[V^] contains 
or higher order terms; the analysis in the real case is similar. We introduce u{r) = u{r) — u~ ^ 
and denote the relevant eigenvalues as ± iA^, where \r^i > 0 , assuming that \i = 0 ( 1 ). 
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When (i = 1, u^^\r = 0) vanishes for j odd and scales as 

^(i) — Q) ^ 0 xp {—Xrp) sin {Xip + (/)) (A. 13) 

for j even, where 0 is a constant. Since H{r = 0) — H{u~) is a quadratic form in u^^\r = 0), the 
higher order correction to Ai due to the oscillatory tail scales as 

M — Mq ^ exp {—2Xrp) sin {2Xip + (^-14) 


where is a constant. Hence c ^ At is a decaying oscillatory function of p; this behavior is known 
as collapsed snaking (CS) 0 . 

For d > M. still exhibits oscillatory decay with p, but with a different scaling from Eq. (jA.14[) 
unless d — 1 ^ 1. The contribution from JC is monotonic in p and scales as 

1C - {d-l)p-\ (A.15) 


If Xr = 0(1) and d — 1 = 0(1), then c ^ Ad JC decays monotonically in p for p 1 . However if 
either A^^lord—l<Cl, then there exists a critical radius p^^^ ^ 1 such that c is monotonic 
in p for p > p^^^ and oscillatory in p for p < p^^^; this is radially collapsed snaking (RCS). The 
critical radius p^^^ satisfies d{A4 + /C)/dp = 0 , i.e. 

exp (— 2 Arp) ^ (d — l)p“^. (A.16) 


In 0an instance of RCS is found numerically at d = 2 and A^ ^ 1 in the 1:1 FCGLE, or specifically 
Eq. (j2.9p . It is an open question whether A^ ^ 1 is possible in Eq. (IV2I1 . In the d — 1 ^ 1 case, 
denoting e = d — 1 , we can solve Eq. (jA.16p asymptotically, obtaining 

= ^(-loge + 21og(-loge)) +0(1). (A.17) 

To take a concrete example, we consider the steady-state axisymmetric quadratic-cubic 
Swift-Hohenberg equation (SH23) in d dimensions, namely 


0 = 




d — 1 
r 



2 


u - 1 - h 2 U^ 


(A.18) 


where 62 and 7 are constants. At 62 = 2.7, the Maxwell point occurs at 7 m = —0.62 between 
u~ = 1.8 and = 0, and the real part of the unstable eigenvalues of u~ is A^ = 0.37. Figure lA.241 
shows the logarithm of the distance from the Maxwell point, I7— 7 m|, as a function of the front radius 
p computed from numerical continuation of axisymmetric fronts between at d — 1 = 0 (blue) and 
d — 1 = 10 “^ (red). In the d — 1 = 0 case, the branch undergoes CS with decay exponent — 2 A^. 
In the d — 1 = 10 “^ case, the branch undergoes RCS with oscillatory and monotonic dependence 
on p respectively for p < p^^^ and p > p^^^ . The predicted critical radius from Eq. (jA.17p is 
pRCS _ 26 + 0(1), while the predicted radius at which 7 = 7m is p = 22 + 0(1), both of which 
agree well with the numerical computations. 
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